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ABSTRACT 

Using a sample of nearly 20,000 massive early-type galaxies selected from the Sloan Digital Sky 
Survey, we study the color-magnitude relation for the most luminous {L > 2.2L*) field galaxies in the 
redshift range 0.1 < z < 0.4 in several colors. The intrinsic dispersion in galaxy colors is quite small 
in all colors studied, but the 40 milli-mag scatter in the bluest colors is a factor of two larger than the 
20 milli-mag measured in the reddest bands. While each of three simple models constructed for the 
star formation history in these systems can satisfy the constraints placed by our measurements, none 
of them produce color distributions matching those observed. Subdividing by environment, we find 
the dispersion for galaxies in clusters to be about 11% smaller than that of more isolated systems. 
Finally, having resolved the red sequence, we study the color dependence of the composite spectra. 
Bluer galaxies on the red sequence are found to have more young stars than red galaxies; the extent of 
this spectral difference is marginally better described by passive evolution of an old stellar population 
than by a model consisting of a recent trace injection of young stars. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: photometry - 
galaxies: statistics - galaxies: fundamental parameters 



1. INTRODUCTION 

The presence of a tight correlation between the 
rest-frame optical color and luminosity of early- 
type galaxies, the so-called color-magnitude rela- 
tion (CMR) or red- sequence, is well esta blished 
in the literature fVisvanathan fc Sa ndagcj 
Larson. Tinslev. fc Caldwelll ll^SOt iL ugE 
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1977 



Zevf. Whitmorc. & Lcvison 1991: Terlevich ct al. 199^ 
Blant on ct al. 2003a; Hogg ct al.. ,2004; Baldrv et all 
2004 ILoDCz-Cruz. Barkhousc. fc Yeel l2004l ). The rela- 



tion implies that the properties of the stellar populations 
in early-type galaxies are strongly dependent on the 
total stellar mass of the system. The tight correlation 
observed along t his r el ations hi p (e.g. ^2Z^. 
^isvanathan fc Sa ndagd Il977t IBower. Lucev. fc EUisI 
11992) requires a strong coherence between both the age 
and metaUicity of stellar populations present in early- 
type galaxies. Measurements of the slope, zero-point, 
and dispersion around the color-magnitude relation 
probe the star formation histories of early-type galaxies 
and can provide physical insight into galaxy formation 
and evolution. 

One possible explanation for the slope of the CMR is 
that massive, and thus more luminous, galaxies retain 
more metals than less massive ones. The deeper po- 
tential wells of massive galaxies may allow few er met- 
als to be expelled in superno v a driven winds dLarson 
lArimotofc Yo^hilllTgSTt IMatteucci fc Torna,mbf 
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[1987t iBr essan. Chiosi fc Tantald Il99(il) . Alternatively, 
iKauffman n fc Charlod l)1998D have shown that the color- 
magnitude relation can be reproduced in hierarchical 
merging models if the epoch of mergers occurs at very 
early times. Since the age of a stellar population is 
degenerate with its metaUicity on the CMR, however, 
metaUicity may not be the sole parameter behind the re- 
lationshin J^^IT97llT97arO'Corme]in 98ntm,osen 98,'t 



. ^^^j^^^^^^- - — — 

IWorthev et al.lll99,'^^ . While both age and metaUicity ef- 
fects may create the CMR, studies of the evolution of the 
relationship toward high-redshift argue that the relation 
is driven primarily by m etaUicity and cannot be ge ner- 
ated by age effects alone l)Kodama fc Arimotolll99'n) . 

The color- magnitude relation of galaxies in clusters has 
been studied extensively. Using high precision U and 
V photometry. Bower. Lucey. fc Ellis ( 1992,) determined 
the dispersion around the CMR for early-type galax- 
ies in Coma and Virgo to be small, 5{U — V) ^40 
milli-mag (mmag). This small scatter implies early- 
type galaxies in clusters must be uniformly old with 
formation epochs earlier than a z ~ 2. 'Ellis et aQ 
lfT99l found S{U - v) 100 mmag for clusters at 
z « 0.54; this scatter, combined with the early age 
of the Universe at the observed epoch, indicates that 
elliptical galaxies must have stopped forming stars in 
these clusters by a z 2.7 for a simple approximation 
to the color evolution of stellar populations. The clus- 
ter MS 1054-03 at a z = .83 has a_29mniag disper- 
sion m theU-B c o lor fm^ Dokku m et al.ll20flfD while 
Ivan Dokkum et all l|2001() and iBlakeslee et all l|2003D 
showed that the scatter in single clusters is only 40 and 
30 mmag at z ^ 1.3. These small dispersions, seem- 
ingly independent of redshift, indicate that the coor- 
dination of star formation between galaxies in clusters 
is quite strong. Using a sample of 158 clusters with 
0.06 < z < 0.34 drawn from the Ear ly Data Release 
of the Sloan Digital Sky Survey rSDSSV IAndreoTil (12001 



2 



Cool et al. 



found that the CMR between clusters is very homoge- 
neous. Locally, the CM R is universal between clusters of 
various masses as well (jMcIntosh et al.ll2005j) . 

While the mean color of red sequence galaxies in over- 
dense regions is only slightly redder than for similar 
galaxies in less dense regions an d the slope of the CMR 
is independent of environment l)Hogg et al J 12004(1 . sev- 
eral studies show that the scatter around the CMR is 
dependent on gala xy environment. A re-analysis of the 
data presented in iSandage &: Vi svanathan (1978) indi- 
cates the dispersion around the CMR in clusters is 2- 
cr sm aller than that for galaxies in g roups or in the 
field l|Larson. Tinslev. fc CaldwelllllQSnli . The scatter in 
galaxy colors in the core of the cluster CL 1358-1-62 at 
z = 0.33 have 5{B — V) ^ 22 mmag while early-types 
at large radii (i? > lAh^^ Mpc) have nearly double 
that value, indicating that the dispersion of galaxy col- 
ors in the field m ay be larger th an that measured in 
clusters (^an Dok kum et al.lll998fl . Early-type galaxies 
in the Hubble Deep Field have a scatter of 120 ± 60 
mm ag in rest-frame U — V with an average redshift of 
0.9 lIKodama. Bower, fc Bel]llT999l) . While this measure- 
ment is rather uncertain, it is larger than that measured 
in clusters, further indicating that field galaxies may have 
less coordinated star formation than galaxies in clusters. 

A majority of the work on the CMR has focused on 
early-type galaxies in clustered environments with lumi- 
nosities near L*. This is natural, as statistically signifi- 
cant samples are more easily selected in these over-dense 
regions and L* galaxies are more common than their 
more massive analogs. In this work, we extend the analy- 
sis of the color-magnitude relation to a sample of nearly 
20,000 very luminous red field galaxies from the Sloan 
Digital Sky Survey. These galaxies are not selected to re- 
side in clustered environments and are chosen to be quite 
luminous with L > 2.2 L*. We create two galaxy sam- 
ples from the SDSS spectroscopic data. At low-redshifts 
(0.1 < z < 0.2) we use galaxies from the entire survey re- 
gion while at higher redshifts (0.3 < z < 0.4), where the 
galaxies are apparently fainter, we select galaxies from a 
270 deg^ region which has been imaged an average of 10 
times and up to 29 times. The repeated imaging of our 
moderate-redshift galaxies makes comparisons between 
the galaxy samples possible with similar fidelity. Our 
sample allows us to probe a new area in parameter space 
— luminous galaxies in the field — and compare our 
results to those found for galaxies in clustered environ- 
ments. 

The paper is organized as follows: §2 describes the 
Sloan Digital Sky Survey and the selection for the galax- 
ies used here. In §3, we present our measurements of the 
slope and scatter of the color-magnitude relations. We 
discuss simple star formation history models in §4 and 
use our measurements of the scatter around the CMR 
as a constraint on the evolutionary history of massive 
early- type galaxies in the field. In §5, we discuss the role 
environment plays on the CMR of galaxies. We construct 
the composite spectra of galaxies across the red sequence 
in §6 before closing in §7. Throughout this work, we use 
Ho = 100 /i km s"^ Mpc-i and (f],„, f^A) = (0.3,0.7) to 
calculate luminosities and distances. When calculating 
ages of stellar populations, we set h = 0.7. All magni- 
tud es used here have bee n corrected for reddening using 
the lSchlegel et al.l l|1998D dust maps. 



2. DATA 

2.1. The Sloan Digital Sky Survey 

The Sloan Digital S ky Survey llYork et alJ l200r ; 
'Stou ghton et alJ |2002b'; 'Abazai ian et all 120031 12004a , 

2005) is imagi ng tt st cradians of the n orthern sky through 
5 passbands fFukugit a et all I1996^ . The imaging is 
conducted with a CCD mosaic in drift-scanning mode 
llCunn et all IT991 on a dedicated 2.5m telescope lo- 
cated at Apache Point Observato ry in New Mexico . 
After the images are processed llLunton et al.l 120011: 
iStoughton et al. 2002; Pier et al. 2003") and calibrated 
((Hogg, et al. -2001: Smith et al.,2002:,Ivezic et al. 2004) . 
targets are selected for spectroscopy with two double- 
spectrographs mounted on the same telescope using 
an automated spectro scopic fiber assignment algorithm 
(|Blanton et alJl2003bD . 

Two spectroscopic galaxy samples are created us- 
ing the SDSS ima ging. The MAIN galaxy sample 
l|Strauss et al.ll20d^ is a complete, flux limited, sample 
of galaxies with r < 17.77. This cut is nearly five magni- 
tudes brighter than the SDSS detection limit of r ~ 22.5, 
and thus the photometric quantities for these galaxies 
are measured with signal-to-noise of a f ew hundred. The 
Luminous Red Galaxy (LRG) sample ((Eisenstein et alJ 
1200 ID selects luminous early-type galaxies out to z ~ 0.5 
with r < 19.5 using several color-magnitude cuts in g, r, 
and i. The combination of these two samples allows us to 
study the broadband colors of luminous field galaxies at 
moderate-redshifts with statistically significant samples. 

While the focal point for the SDSS is a contiguous 
survey of the Northern Galactic Cap, the SDSS also con- 
ducts a deep imaging survey, the SDSS Southern Survey, 
by repeatedly imaging an area on the celestial equator 
in the Southern Galactic Cap. Currently, objects in the 
270 deg^ region have been imaged an average of 10 times 
and up to 29 times, resulting in improved photometric 
precision for faint galaxies. Objects detected in each 
observational epoch were matched using a tolerance of 
0.5 arcseconds to create the final coadded catalog used 
here. The photometric measurements from each epoch 
were combined by converting the reported asinh magni- 
tudes (Lupton et al. 1999) into fiux and then calculating 
the mean value. Errors on each parameter are simply 
the standard deviation of the fiux measurements. The 
Southern Survey is an ideal region to compare the prop- 
erties of faint galaxies at moderate-redshift with brighter 
low-redshift analogs drawn from the entire SDSS survey 
area with photometry of similar fidelity. 

2.2. Galaxy Photometric Properties 

Several methods are used to measure galaxy fluxes in 
SDSS. Below, we briefly describe the two flux measure- 
ments used throughout this paper. The Petrosian ratio, 
TZp, the ratio of the local surface brightness at radius 6 
to the average surface brightness at that radius, is given 

by 

T^p(O) = -, , (1) 

Jo dr'2^e'i{e') I [^02] 

where I{9) is the azimuthally averaged surface bright- 
ness profile of a galaxy and a\o and ahi are chosen to 
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Fig. 1. — Sample definition for the galaxies used in tliis paper. Tlie left panel shows the low redshift selection regions; contours encapsulate 
25%, 50%, 75%, and 95% of the considered galaxies. The right panels shows the target selection for our liigh-redshift galaxies overlayed on the LRG 
galaxies in the SDSS S outhern Survey . In each panel, the solid line shows the selection cut which approximates a cut at constant stellar mass using 
the relations of ,,Bell de Jonel 1200011 . The regions defined by the dashed lines were used to create control samples which do not have a strong cut 
across the red sequence to test the effects of sample selection on our scatter measurements. The arrow in the lower right of the left panel shows the 
reddening vector. Note the dearth of very luminous blue galaxies which could contaminate our red-sequence samples if heavily reddened. 



be 0.85 and 1.25 for SDSS. The Petrosian flux is given 
by the flux within a circular aperture of 29 p, where 9p 
is the radius at which TZp falls below 0.2. In SDSS, Op 
is determined in the r band then subsequently used in 
each of the other bands. This flux measurement con- 
tains a constant fraction of a galaxy's light, independent 
of its size or distance, in the absence of seeing. More 
details of SDSS Pet rosian magnitudes can be found in 
IBlanton et all l|2001|) . and lStrauss et all l|2002D . 

For each galaxy imaged by SDSS, two seeing-convolved 
models, a pure dc Vaucoulcurs (1948) profile and a pure 
exponential profile, are fit to the galaxy image. The best- 
fit model in the r band is used to measure the flux of a 
galaxy in each of the other bands. These model magni- 
tudes are unbiased in the absence of color gradients and 
provide higher signal-to-noise ratio colors than Petrosian 
colors. A more complete descrip t ion of model magni- 
tudes is given in IStouehton et all lj2002b|) . Throughout 
this paper, we use Petrosian magnitudes when calculat- 
ing galaxy luminosities while model magnitudes are used 
to determine galaxy colors. 

2.3. The k- Correction 

We cal culate photometric fc -corrections using the 
method of lBlanton et alJ l)2003cj) (kcorrect v3_2). This 
program constructs a linear combination of carefully 
chosen spectral templates in order to best match the 
observed photometry at the measured redshift of the 
galaxy. The rest-frame colors of these best-fit spectra are 
used to derive corrections to the observed galaxy colors. 
With this method, it becomes simple to transform be- 
tween bandpasses and, if necessary, photometric systems 
given a detailed understanding of the filter characteris- 
tics. After calculating fc-corrections with kcorrect, we 
remove any mean A:-corrected color versus redshift trends 
with a low order polynomial to account for the passive 
evolution of stellar populations between observed epochs. 
This evolutionary correction is normalized to a redshift 
of 0.3, near the median of the LRG redshift distribution. 
The color evolution correction is small (less than 2% for 



TABLE 1 

Effective 
Wavelengths 
OF Shifted 
Bandpasses 



Filter 




(1) 




O.Wg 


4026 


0.16^ 


5307 


0.16j 


6441 


o.ie^ 


7687 


0.37^ 


3408 


0.37^ 


4494 


0.37j 


5454 


0.37^ 


6509 



all bands in both samples) and does not affect any of the 
results presented here. We further assume a galaxy dims 
by one magnitude per Az = 1 change in redshift in the 
g-band due to the passive evolution of its stellar popu- 
lations. Again, this correction is normalized to z = 0.3. 
In the remainder of this paper, all of the color and lumi- 
nosity measurements refer to quantities which have been 
adjusted for the fc-correction and evolution. 

Throughout this paper, we work in the '^'^^g, "-^^r, 
o.i6j^ 0.16^ a,nd °-^'^g, ^ ^^r, "■^'^i, "-^^^ systems for the 
low- and moderate-redshift galaxy samples, respectively. 
Here, the notation indicates the reported AB mag- 
nitudes are derived through a standard g filter that has 
been blueshifted by zq. For a filter system which has 
been shifted by zq, the fc-correction for galaxies at that 
redshift is trivial, — 2.51og(l + zq), and is independent 
of the galaxy spectral energy distribution. In choosing 
shifts which nearly match the median redshift of each 
sample, we minimi ze the uncertainty intr oduced through 
the fc-corrections ([Blanton et alJl2003c(l . Table 1 lists 
the effective wavelengths for each of the bandpasses used 
here. 
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2.4. Sample Construction 

In order to investigate the color-magnitude relation of 
massive galaxies in the range z ~ 0.1 — 0.4, we create 
two samples from the SDSS spectroscopic database to be 
analyzed separately; Figure Q shows the selection cuts 
used to define each of these samples. From the MAIN 
sample of galaxies, we extracted 16622 galaxies in the 
redshift range 0.1 < z < 0.2 satisfying the rest-frame 
color cut 

°-^^{g~i) > 8.374 + 0.32 (Mg- 5 log /i) (2) 

which approxiniates a cut of constant stellar mass 
ijBell fc de Jond l2000j) . This selection criterion results 
in a cut diagonally across the red sequence, which could 
bias our measurements of the scatter in red sequence col- 
ors and will greatly affect any slope measurement of the 
color-magnitude relation. In order to avoid biased slope 
measurements and test the consistency of our scatter 
measurements, we also construct a sample of 5100 galax- 
ies in the luminosity range —22.5 < Mg — 5 \ogh < —21.5 
with '^'^^{g — i) > 1.46 for comparison (shown by the 
dashed region in Figure P). 

We further select 2782 LRG galaxies from the SDSS 
Southern Survey in the redshift range 0.3 < z < 0.4 
which satisfy 

°-3^(r - z) > 7.769 + 0.32 (Afg - 5 log /i) (3) 

and have been imaged 6 or more times. Equation 3 is 
the same rest-frame color cut applied to the low-redshift 
galaxies (Equation 2) shifted to the '^■^'^(j — z) color — 
the effective wavelengths of the filters that define each 
of these colors are similar, and thus these colors probe 
similar regions in the rest-frame galaxy spectrum (see 
Table 1). Of these 2782 sample galaxies, 62% were 
imaged more than 9 times while 20% of the galaxies 
were imaged more than 12 times in the SDSS South- 
ern Survey. As in the low-redshift sample, we adopt 
a simple constant color cut, 0.37 j-^, _ ^-j 0.825 and 
—22.5 < Mg — 5\ogh < —21.5, to define a subsample 
of 1380 high-redshift galaxies to test the consistency of 
our scatter measurements and allow us to determine the 
slope of the color-magnitude relation without a strong 
cut across the red sequence. 

We find no differences within our quoted error between 
the dispersion around the CMR measured from galaxies 
defined using the stellar mass cut and the constant color 
cut; throughout the remainder of this paper, scatter mea- 
surements will be based upon galaxies drawn in the first 
manner while all slope measurements are referenced to 
samples created using the latter method. FigureElshows 
the redshift distribution of the galaxies used to measure 
the scatter around the color-magnitude relation in this 
paper. 

In detail, the red sequence is populated by elliptical 
and SO galaxies as well as late-type galaxies reddened by 
dust. The arrow in the lower ri ght of Figure H] sh ows 
the reddening vector given by the lO'Donnelil 1119941 ) ex- 
tinction curve. In order for a dusty late-type galaxy to 
scatter into our sample, it must have a large unobscured 
luminosity as we only consider the most luminous red- 
sequence galaxies in this work. There is a clear paucity 
of very luminous blue galaxies in Figure ^ and thus the 
contamination from dusty spiral galaxies is likely small 
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Fig. 2. — Redshift distribution of the galaxies selected for this 
work is shown. The low-redshift sample was drawn from the MAIN 
SDSS galaxy sample while the moderate-redshift sample is composed 
of LRG galaxies selected from the SDSS Southern Survey. The low- 
redshift sample has robust photometry in one observation epoch as 
the galaxies are upwards of five magnitudes brighter than the SDSS 
detection limit while the higher redshift galaxies, which extended two 
magnitudes fainter, have been imaged repeatedly allowing for more 
precise photometric measurements. 



in our sample. Furthermore, lEisenstein et alJ l|2003D con- 
structed the average spectrum of massive galaxies, such 
as those used in our analysis, and found that only 5% of 
the galaxy spectra in the —22.5 < Mr < —22 luminosity 
range were contaminated by emission lines uncharacteris- 
tic of early-type galaxies. In the remainder of this paper, 
we will use the red-sequence and early-type classifications 
synonymously. 

Throughout this paper, we will refer to the galaxies 
considered here as field galaxies as they are not specif- 
ically selected to reside in clustered environments. It 
is important to note, however, that galaxies which are 
of similar color and luminosity as th ose in our sampl e 
tend to lie in o ver-d ense environments ijHogg et alJ2003j) . 
iZehavi et al.l l)2005f ) found the clustering strengths and 
mean separations of LRGs are comparable to those of 
poor clusters or rich groups. Fo r comparison, of the 
LRGs in the volume surveyed by iBahcall et al.l l|2003f) 
to find galaxy clusters with 13 or more detected clus- 
ter members, 16% are projected less than 500 kpc from 
the inferred location of a cluster and within 0.05 of the 
photometric redshift for that cluster. A detailed exam- 
ination of the density dependence of the scatter around 
the CMR is conducted in §5. 
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Fig. 3. — Measured slope of the CMR as a function of the rest-frame 
bandpass. The endpoints of the horizontal error bars mark the effcetive 
wavelengths of the filters used to define the bandpass while the vertieal 
error bars show the measured uneertainty on each slope measurement. 
The spectrum of an early-type galaxy is plotted as a reference. The 
two symbols represents the two galaxy samples - the hollow squares 
mark measurements on the moderate-redshift (0.3 < z < 0.4) sample 
while the filled circles show the data from low-redshift (0.1 < z < 0.2) 
galaxies. 



3. THE COLOR-MAGNITUDE RELATION 

When determining the slope of the color-magnitude 
relation, we only consider the galaxies selected with our 
simple luminosity and constant color cuts as any slope 
measurements based on samples with a strong cut across 
the red sequence will be quite biased. The galaxy sam- 
ples considered in this work only span a small range of 
luminosities and thus are not ideally suited for measure- 
ments of the slope of the CMR - we simply use our 
measurement of this quantity to remove the mean re- 
lationship before determining the scatter around it. The 
slopes and zero-points for the CMR for several colors are 
listed in Table 2. Figure O shows the measured slopes 
of the CMR as a function of the effective wavelengths of 
the two bandpasses used to define the color. The slope 
in bluer colors is more pronounced than that for red- 
der bands which tend to show small or negligible slopes. 
This trend is consistent with a metallicity sequence as 
the primary driver along the red-sequence. As noted 
bv iGladders et all l)1998j) using models from iKodamal 
the flattening of the color-magnitude relation to- 
ward redder colors can also be reproduced by an age 
trend along the red sequence but only if a very spe- 
cific relationship between galaxy mass and star formation 
time holds with very little scatter. In the bluest colors, 
which have been used historically in the literature, the 
slopes measured here are in general agreement with past 
work; for comparison, Hogg et al. ( 2004) report a slope 
of —0.022 mag mag^^ in the — r) color for a large 
sample of galaxies at z ~ 0.1. 

Figure ^illustrates the typical scatter about the mean 



color-magnitude relation for the "■^^(r — i) and '^■^^(j. _ 
colors; the measured dispersion about this relationship 
is quite small. The second column of Table 3 lists the 
measured dispersion around the CMR for each of the 
colors used in this study. For each color studied, we 
fit a Gaussian to the observed color distribution using 
Poisson errors on each color bin and define the disper- 
sion of the best-fit Gaussian to be the scatter in galaxy 
colors on the red sequence. By fitting the histograms 
rather than calculating the standard deviation directly, 
we place more weight on the core of the distribution com- 
pared to the wings. Thus, our scatter measurements will 
not be heavily affected by any non-red-sequence inter- 
lopers in the wings of the color distributions. As this 
method could potentially depend on the binning used, 
we are careful to vary the bin size for each fit and ensure 
the measured scatter is not dependent on our choice of 
bin size. It should be noted that past work has used the 
standard deviation to quantify the scatter in galaxy col- 
ors. While the standard deviation and Gaussian widths 
may give systematically different dispersions in the case 
of non-Gaussian distributions, Figure 0] shows the color 
distribution on the red sequence is well described by a 
Gaussian and thus the two measurements a re compara- 
ble. Bootstrap resampling of the data (e.g. lPressll200^ 
was employed to determine the uncertainty on our scatter 
measurements. It is important to note that the scatter 
measurements in the second column of Table 3 represent 
upper limits to the intrinsic scatter in galaxy colors as 
these values are not corrected for systematic or experi- 
mental effects, which could be significant. 

In order to investigate the level of skewness of the ob- 
served color distribution, we define /3„ to be the absolute 
value of the color difference between the n% quantile and 
the median of the distribution. The level of asymmetry 
in the distribution is then estimated using 
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/3io - A 



90 



/3io + f3'. 



'90 



(4) 



For this statistic, a maximally skewed distribution would 
have 7 = ±1 such that a tail of galaxies to the red would 
result in 7 < 0. A symmetric distribution would have 
7 = 0. We adopt this approach rather than a traditional 
third moment calculation as the calculation of higher 
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Fig. 4. — Scatter around the color-magnitude relation for both 



galaxy samples. The 



i) and 



residual colors, normal- 



ized to the characteristic color of a galaxy with Mg — 5 log /i — —21.8 
arc shown for the low and intermediate redshift galaxies respectively. 
Note that the dispersion of this relation is quite small. The dashed line 
is the best-fit Gaussian to the data. Table 3 lists the measured scatter 
around the CMR for all of the colors measured in this paper. 
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TABLE 2 

Slope of the Color-Magnitude Relation 







Total Sample 


Field Sample 


Cluster Sample 


Color 




Zero-point^ 


olope 


Zero-point^ 


olope 


Zero-point^ 


olope 






mag 


mmag/mag 


mag 


mmag/mag 


mag 


mmag/mag 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


0.16(^_ 


r) 


1.162 


-15.6 ± 5.2 


1.162 


-10.5 ±6.7 


1.164 


-12.4 ±7.8 


0.16(^ _ 


i) 


0.449 


-8.10 ±2.6 


0.449 


-6.7 ±3.4 


0.451 


-7.6 ±4.3 


0.16(j_ 


^) 


0.329 


3.50 ±6.4 


0.329 


4.2 ±8.2 


0.329 


5.5 ± 10.4 


0.16(3 _ 


i) 


1.611 


-26.1 ± 5.5 


1.610 


-21.5 ± 7.4 


1.612 


-22.9 ±8.2 


0.37(3 _ 


r) 


1.750 


-23.3 ± 12.1 


1.750 


-15.5 ± 16.1 


1.754 


-10.6 ± 19.5 


0.37(j. _ 


i) 


0.613 


-3.8 ±4.3 


0.612 


0.41 ±6.2 


0.614 


-3.6 ±8.4 


0.37 (j _ 


2) 


0.379 


5.0 ±6.1 


0.377 


0.58 ± 7.7 


0.381 


7.21 ±9.9 


0.37(j. _ 




0.991 


1.8 ±8.5 


0.980 


3.4 ±3.5 


0.983 


8.1 ± 1.4 



^The zero-point is defined at Mg — 5 log/i = —21.8. 



moments of a distribution is strongly dependent on the 
wings of the distribution where contamination may be 
important. For the histograms shown in Figure 2) we 
find that 7 = -0.027 ±0.019 for the ^ ^^(r - color and 
7 = -0.064 ± 0.048 for the °-3^(r - i) color. Both of 
these measurements reflect the slight over-abundance of 
red galaxies compared to a Gaussian in both panels of 
Figure El This red tail is due to the diagonal cut across 
the red sequence imposed by our stellar mass selction 
criterion; red galaxies are selected to lower luminosi- 
ties than blue galaxies. Since the luminosity function of 
galaxies is steeply rising toward lower luminosities in this 
range, more red galaxies are selected than blue galaxies. 
For comparison, we measure 7 — 0.001 ± 0.029 for the 
°-i^(r-i) color and 7 = -0.005 ±0.057 for the "■^^(r-i) 
color when the simple luminosity cuts are used to select 
our sample galaxies. 

Since our intermediate redshift sample is drawn from 
regions that have been repeatedly imaged by the Sloan 
Digital Sky Survey, we can assess the level to which the 
dispersion we observe in the color-magnitude relation is 
dependent on measurement noise. If the scatter is con- 
taminated by measurement noise, we would expect the 
squared scatter to decrease inversely with the number of 
observations : 

on A 



obs 



intrinsic 



(5) 



Figure shows the results of this test; the methods used 
for the ensemble sample were used in the same manner 
to calculate the scatter and error for each bin. In each 
bin of A^-measurements, we include galaxies which were 
imaged a total of N times, and thus each bin contains a 
unique subsample of galaxies. For the moderate-redshift 
sample, the third and fourth columns in Table 3 lists 
the instrumental contribution and intrinsic galaxy scat- 
ter predicted by extrapolating the fits to Equation 2 to 
infinite observations for each of the three colors. 

While the intrinsic dispersion around the color- 
magnitude relation at moderate redshift can be extracted 
from the measured scatter based on subsamples of galax- 
ies with different number of measurements, the low- 
redshift galaxies do not allow for this test as most galax- 
ies in this sample have been imaged only once. In order 
to estimate the intrinsic scatter in galaxy colors for this 
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Fig. 5. — Relationship between the number of observations and the 
measured scatter around the color-magnitude relation. The dotted line 
shows the weighted mean of the data while the dot-dashed line marks 
the best fit line to the data. The data for each color are not statistically 
consistent with a constant value but are well described by a line, as is 
expected. We extrapolate the best fit line to infinite measurements to 
find the intrinsic scatter in galaxy colors; these values are listed in the 
fourth column of Table 3. 



sample, we identified 5839 of our low-redshift galaxies 
that have been imaged more than once under photomet- 
ric conditions. We adopt the average root-mean squared 
variation of each photometric measurement as the in- 
strumental error for each color of interest. We find the 
mean measurement uncertainty to be 17.6, 15.1, and 20.9 
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TABLE 3 

Scatter Around the Color-Magnitude Relation 



Color 




Measured Scatter 


Instrumental Scatter 


Intrinsic Scatter 


Field^ 


Clustered^ 






mmag 


mmag 


mmag 


mmag 


mmag 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


0.16(g_ 


r) 


39.6 ±0.4 


17.6 ±3.6 


35.4 ±3.7 


40.5 ±0.4 


35.9 ±0.8 


0.16(-^ _ 


i) 


18.7 ±0.2 


15.1 ±3.9 


11.0 ±4.0 


19.1 ± 0.2 


17.3 ±0.3 


0.16(j_ 


z) 


26.6 ±0.4 


20.9 ±4.2 


20.9 ±4.2 


26.7 ±0.5 


23.4 ±0.8 


0.16(g_ 


i) 


40.0 ±0.4 


21.2 ±4.1 


34.4 ±4.1 


42.0 ±0.5 


36.1 ±0.7 


0.37(3 _ 


r) 


48.0 ±1.2 


20.4 ±6.1 


43.5 ±6.2 


48.8 ± 1.4 


47.3 ±2.2 


0.37(-^ _ 


i) 


19.7 ±0.5 


12.3 ±3.1 


14.6 ±3.1 


20.1 ± 0.6 


18.9 ±1.0 


0.37 (i_ 


z) 


27.9 ±0.7 


21.6 ±4.5 


17.6 ±4.6 


29.2 ± 0.8 


24.4 ±1.4 


0.37(j, _ 


z) 


35.9 ± 1.1 


18.4 ±5.2 


30.8 ±5.2 


38.1 ± 1.2 


28.4 ±2.1 



"Field and clustered values are measured scatter, not the intrinsic scatter in the galaxy colors. 



mmag in the '^'^^{g — r), '^■^^(^r — i), and — z) bands, 
respectively. The internal dispersion in each color is the 
difference, in quadrature, between the measured scatter 
and inferred mean uncertainty in our single pass photom- 
etry. For the low-redshift galaxies, the third and fourth 
columns of Table 3 lists the observational error and in- 
trinsic scatter determined by this correction. While our 
estimate of the instrumental contribution to the observed 
color dispersion corrects for uncorrelated errors on each 
independent measurement, we do not correct for system- 
atic errors common to all observations which could cause 
an overestimation of the intrinsic dispersion in galaxy 
colors. 
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Fig. 6. — Intrinsic scatter in the color-magnitude relation as a 
function of the observed bandpass. The error bars and symbols are 
defined as in Figure|3] The spectrum of an early-type galaxy is shown 
for comparison. The scatter in redder bands tends to be smaller than 
in blue bands. This is not unexpected, as small changes in a galaxy's 
age, metallicity, or dust attenuation affect the blue portion of a its 
spectrum more than the red. 

Figure ini shows our scatter measurements as a function 
of the effective wavelength of the filters used to construct 



each color. The scatter about the OMR is clearly wave- 
length dependent; the bluest colors, which are most af- 
fected by small changes in the age, metallicity, or dust 
attenuation of the system, have systematically larger dis- 
persions than the reddest colors where changes in metal- 
licity and age have less impact. 

4. STAR FORMATION HISTORY OF MASSIVE GALAXIES 

Here, we consider three simple star formation histo- 
ries and derive the range of permitted parameters for 
these histories following a method similar to that of 
Ivan Dokkum et all ljl998 ^. In these toy models, we will 
assume that only age variations between red-sequence 
galaxies creates the scatter in early-type galaxy colors. 
In reality, the scatter across the CMR is likely driven 
by some combination of age, metallicity, dust, and pos- 
sibly a-enhancement variations across the red-sequence. 
These models will thus illustrate the broadest range of 
ages allowed by our measurements; the inclusion of any 
variation in the other galaxy properties will only con- 
strain the age more strongly. 

In this section, we will primarily consider the high- 
redshift sample, as these are the intrinsically youngest 
galaxies in our sample and hence should provide the 
tightest constraints on the past star formation history 
of massive early- type galaxies. An identical analysis of 
the low-redshift sample provides results consistent with 
those presented below. In all three model cases, the low- 
redshift galaxies provide looser constraints than the high- 
redshift sample; star formation is allowed to continue to 
slightly later times. 

First, we consider a star formation history composed of 
a single delta-function burst. The probability of a burst 
occurring is distributed uniformly between the earliest 
allowed time for galaxy formation, tg, set to z = 10, and 
^maxj the time at which the youngest early- type galaxies 
form. Using a grid of Monte-Carlo simulations with 5000 
galaxies per realization, we create sam ples of galaxy spec- 
tra sy nthesized using the method of Bruzual fc C hariot! 
lj2003j) with solar metallicity, a Salaeter (1955,) initial 
mass function, and varying values of fmax- We find 
that realizations allowing galaxies younger than 8.8 Gyr 
(z ^ 1.4) generate color distributions with standard de- 
viations larger than that measured. 

We perform a similar test using a scheme in which 80% 
of the stars in a galaxy are formed in a single burst at 
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z = 10 and 20% are formed in a second burst which 
occurs at a random tir ne between z = 10 and a min- 
imum age, as done by Ivan Dokkum et all l(T99^ . For 
this scenario, the secondary bursts must occur between 
10 > z > 0.9 to reconstruct the observed scatter in the 
color-magnitude relation. For bursts that generate less 
than 5% of a galaxy's stars, the time scale for the sec- 
ondary burst is unconstrained - small amounts of star 
formation can occur at late epochs without violating the 
constraints placed by our measurements. 

Finally, we create a grid of models in which star for- 
mation begins at z = 10 and continues at a constant rate 
until it is truncated at tf. For each realization, we allow 
t{ to range from z = 10 to a critical cutoff time at which 
all galaxies have stopped forming stars. Star formation 
in this manner must have ended by z ~ 0.8 in order to 
reproduce the observed scatter in the color-magnitude 
relation at z = 0.37. 

The allowed range of parameters and the distribution 
of galaxy colors found in each model are shown graph- 
ically in Figure [7| To allow comparison with Figure 01 
typical observational errors have been included in the 
color distributions in the lower panels of the figure. Inter- 
estingly, all of these models predict similar scatter trends 
to those seen in Figure|B| In all three cases, the scatter in 
the — r) color is nearly double that predicted for the 
0.37 J- J, _ j-j g^j^^ "-37 (^j _ colors. While each of the three 
models succeeds in reproducing the observed scatter in 
early-type galaxy colors, none of the models adequately 
recreates the color distribution measured from real galax- 
ies, a sign that the models chosen do not perfectly track 
the true evolution of the star formation history of early- 
type galaxies. The 7 parameter, as defined in Equation 
4, for each of the model distributions is shown in Fig- 
ureEl The three models produce color distributions with 
7 (0.165, 0.508, 0.313) compared to the value of -0.064 
calculated for the observed sample of galaxies at the same 
redshift. The difference between the observed skewness 
and that determined from our simple models reiterates 
the mismatch between the observed and predicted color 
distributions. In our simulations, we assume all galaxies 
have the same metallicity and no dust attenuation. Also, 
the onset of star formation is uniform in the second two 
models. In reality, these parameters are not constant, 
and thus the sharp truncations seen on the red edges 
of model color distributions will be blurred by the addi- 
tion of other complications to the models. However, any 
of these corrections will increase the observed scatter in 
galaxy colors thus driving the last epoch of star forma- 
tion needed to satisfy our constraints to earlier times. 

5. COMPARISON WITH CLUSTERED ENVIRONMENTS 

A majority of the work on the color- 
magnitude relation of early-type galaxies has fo- 
cused on cluster galaxies at various rcdshifts. 
While som e work ( e.g. ba ndage & Visvanathan 
19781 ILarso n. Tinsl eY, fc Caldwell! I198nt 

Kodama. Bower, fc Bell 1993) tackles the problem 
in the field, these studies suffer from small numbers of 
galaxies. The measurement of the scatter around the 
CMR presented here, on the other hand, focuses on the 
most massive field galaxies at moderate rcdshifts using 
a statistically significant sample. We can now, for the 
first time, statistically compare the dispersion around 



the CMR in both clustered and field environments. 

In order to explore the environmental dependence of 
the color-magnitude relation within our data set, we use 
SDSS imaging of 5305 deg^ to identify a comparison 
sample of 16 million normal galaxies down to r = 21. 
At the location of each LRG galaxy in our sample, we 
count the number of neighboring galaxies within a lh~^ 
Mpc (proper) aperture which have colors similar to those 
expected for an L* galaxy on the red sequence, and 
luminosities in the M* - 0.5 to M* + 0.4 range; the 
details of the luminos ity and color cuts are given in 
lEisenstein et all l|2005j) . 

Using the number of neighbor galaxies as a proxy for 
the environment of each LRG in our sample, we di- 
vide each of our galaxy sets into a clustered subsample 
and a low density subsample. We define LRG galaxies 
with more than 5 (7) neighboring galaxies in the low- 
(moderate-) redshift samples to be in clustered environ- 
ments while a field subsample is comprised of galaxies 
with fewer neighbors. This threshold was chosen such 
that the clustered sample in both redshift bins composes 
30% of the total galaxy sample at that redshift. To place 
these values in physical context, at z = 0.16 (z — 0.37) 
there are 25.3 (170) tracer galaxies per square degree. 
Given the angular diameter distance of 398 (712) 
Mpc to that redshift, this yields a total of 1.65 (3.46) 
background galaxies per aperture on average. The aver- 
age number of neighbors in the field sample is 2.3 (3.9) 
galaxies while the clustered galaxies have a mean of 9.1 
(11.4) neighbors. Thus our clustered samples are over- 
d ense by a fac tor of 4 (3) compared to the field galaxies. 

iHogg et al.l |2004) found that the slope and zeropoint 
of the color-magnitude relation do not depend strongly 
on the local galaxy environment. We find the same re- 
sult; Table 2 lists the measurements of the slope and 
zeropoint of the relationship for the ensemble sample as 
well as for the field and clustered subsamples. The vari- 
ation between the two samples is small. In the following 
calculations, we use the ensemble slope and zeropoint 
measurements as these are measured with the highest 
signal-to-noise. The results are unaffected if the individ- 
ual measurements of the slope and zeropoint are used for 
each subsample. 

The fifth and sixth columns of Table 3 list the mea- 
sured scatter in galaxy colors for the field and clustered 
galaxies respectively. In all of the bands, the scatter 
for the field subsample is larger than that in clustered 
environments. If we assume the scatter introduced by 
instrumental effects is the same for both the field and 
clustered subsamples, the same correction determined for 
the ensemble galaxy sample can be applied to these sub- 
samples, as well. With this correction, we find that the 
scatter in galaxy colors in dense environments is 11 ± 2% 
smaller than for that of isolated systems in both redshift 
intervals studied here. 

Our clustered galaxies reside in a range of densities, not 
simply rich clusters, as have often been used for similar 
studies, and thus we likely underestimate the significance 
of the difference in scatter between the field and strongly 
clustered environments. Furthermore, our field sample is 
not completely composed of isolated galaxies as the aver- 
age number of neighbor galaxies in this sample is larger 
than the average density of background galaxies. Again, 
this leads to an underestimation of the true difference in 
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Fig. 7. — Reproduction of the observed scatter in rest-frame color of carly-typc galaxies using Monte-Carlo simulations following a method 
similar to that of^an Dokkum et al. 1 199^. For each simulation, the top panel shows several examples of star formation histories drawn from the 
simulation. The shaded region denotes the allowed epoch of star formation for that simulation. Each of the bottom panels shows the distribution of 
residual colors around the color-magnitude relation for the realization that most closely matches the observed scatter. The 7 parameter, as defined 
in the text, is listed for each of the distributions. For comparison with the color distributions in Figure |^ we have added typical observational 
errors to the color distributions. The dashed line in each panel shows a Gaussian with the same mean and dispersion as the simulated colors. The 
left-most panel illustrates the allowed epoch of star formation for a SFH composed of a single delta function burst. For galaxy formation beginning 
at a redshift of 10, the youngest early-type galaxies must have formed by z ~ 1.13. The middle panel illustrates the allowed range of times for 
a secondary burst which contributes 20% of the stars to the galaxy. This secondary burst must have occurred before a 2: ~ 0.84. Finally, for a 
constant star formation history, all star formation must have ended by z ~ 0.62 to reproduce the observed scatter in the color-magnitude diagram 
at 2 — 0.37. For reference, the lookback time at 2 — 0.37 is 3.89 Gyr. A comparison with Figure |4] shows that the color distributions predicted 
by these models do not match the observed colors of massive galaxies, a sign that the models used here are overly simple. This is echoed by the 
mismatch between the 7 value measure from the data and those measured for each of the model distributions. 



the dispersion around the color-magnitude relation for 
field and clustered galaxies. 

The larger scatter in field galaxy colors supports pre- 
vious claims of an environmental dependence for the 
scatter around the color-mag nitude relatio n for carly- 
type galaxies (e.g. iLarson. Tinslev. fc Caldwelll I198Q[ 
iKodama. Bower, fc Belli Il999() . though the effect mea- 
sured here is weak. Comparisons of the scatter in galaxy 
colors for clustered galaxies often draw on galaxies from 



a single cluster whereas our clustered sample is created 
by combining galaxies from a number of clustered envi- 
ronments. It is possible that the star formation histories 
of galaxies in a single cluster are more highly coordi- 
nated than the star formation histories between different 
clusters, thus the scatter measured using galaxies from a 
single cluster would be smaller than when several clus- 
ters are combined. As the fraction of blue galaxies is 
larger in the field compared to clustered environments, 
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the field subsample may have a larger contamination 
fraction compared the the clustered subsample. This 
would also cause an increased scatter in field galaxies 
colors compared to the galaxies in more dense environ- 
ments. 

6. AVERAGE GALAXY SPECTRUM ACROSS THE RED 
SEQUENCE 

Having measured the scatter in galaxy colors, we next 
investigate the spectral differences between galaxies by 
examining the average spectrum of massive early-type 
galaxies across the red sequence. The composite early- 
type galaxy spectrum has been investigated as a func- 
tion of environmen t, luminosity, redshift, size, and ve- 
locity dispersions ( Eise nstein et al. 2003'; iBernardi et al.l 
1200,'^ : here, we construct the composite spectrum for 
red sequence galaxies of different mean colors. Since the 
galaxies in our low-redshift sample are brighter, and thus 
have higher quality spectra, and are more numerous than 
the higher redshift galaxies, we only consider these galax- 
ies in our analysis. 

To construct the average spectrum of the galaxies in 
each color bin, we first shift the spectra according to the 
observed redshift and then normalize the observed SDSS 
spectrum by the flux at rest-frame 5000 A. We augment 
the SDSS pixel mask to include all pixels within 10 A 
of the bright sky lines at 5577, 5890, 5896, 6300, 6364, 
7245, and 7283 A. Since the galaxies of interest are in the 
redshift range 0.1 < z < 0.2, the masked pixels occur at 
different rest-frame wavelengths for each galaxy, leaving 
the mean spectrum unaffected. We then find the number- 
weighted mean of the spectra, ignoring all masked pixels 
in the mean. 

Galaxies are divided into two distinct color bins - red 
galaxies having ^'^^{g — r) colors between 1-cr and 3-cr 
redward of the mean and a blue sample in the same 
range blueward of the mean. The red galaxy bin con- 
tains 1321 galaxies with a mean "'^^(5 — r) color of 1.20 
and Mg — 5\ogh = —21.3 while the blue sample con- 
tains 1161 galaxies with mean color and luminosity of 
1.09 and Mg - 5 log /i = -21.5. Figure IHIa) shows the 
coadded spectra of the red and blue galaxies. The gross 
properties of each composite spectrum are similar. Fig- 
ure IHl^b) shows the difference between the blue and red 
spectra. This difference spectrum shows clear signs of 
enhanced Balmer absorption in the blue galaxies com- 
pared to the red composite, indicating that, on average, 
the blue galaxies cont ain more young stars, in agreement 
with past studies fe.g. lTraeer et alJl200Ct iBernardi et al.l 
12005(1 . While the red composite shows stronger [N II] 
emission lines relative to the blue composite, no other 
prominent emission lines, such as [Oil], [Olll], and Ha 
are present; this emission line difference is likely due to 
a larger fraction of AGN in the red galaxies compared to 
the blue. 

We consider two possible evolutionary scenarios to ex- 
plain the spectral differences between the blue and red 
galaxies on the red sequence: pure age (i.e. passive) 
evolution of a single age stellar population and a recent 
epoch of star formation mixed with an existing old pop- 
ulation. In the first scenario, we fit each composite spec- 
trum with a solar metallicity, single age, stellar synthesis 
model created using the method of ..Bruzual fc Gharloti 
l|2003(l . The red and blue spectra are best fit by 6.3 Gyr 



and 3.7 Gyr populations, respectively. In the second sce- 
nario, we find the blue composite spectrum is well fit 
by a 6.3 Gyr old stellar population mixed with a 1 Gyr 
population containing 3.5% of the stellar mass and we 
model the red composite with a simple 6.3 Gyr popula- 
tion. The difference spectra for both scenarios are shown 
in Figure IHIb). While each of the two model hypothe- 
ses would suggest quite different evolutionary histories 
for massive early-type galaxies, the predicted difference 
spectra are remarkably similar. The model consisting of 
a recent burst of star formation slightly over-estimates 
the continuum in the difference spectrum near 3900 A 
and thus provides a marginally poorer fit to the data 
than the passive evolution model. 

The residuals between the data and the model differ- 
ence spectra are shown in Figure |H{c). Both of the mod- 
els fail to match the strengths of the CN2 and Mg b 
features near 4140 and 5160 A. lEisenstein eraP pOOl 
interpret these differences between observations and syn- 
thetic spectra as a enhanced a-to-iron ratio present in 
the galaxies but not included in the synthetic models. 
Also, several weak iron features are evident in the two 
model spectra near 4300 and 5300 A which are miss- 
ing in the observed difference spectrum. In a pure-age 
model for the evolution across the red sequence, the en- 
hanced blue continuum provided by the hotter stars in 
the young population dilutes the metal lines seen in the 
older population. This dilution leads to differences in 
the strengths of the metal lines of the two population 
models, thus causing the residual iron features in both 
model difference spectra. The absence of these residu- 
als in our observed difference spectrum indicates that a 
pure-age model for the spectral differences between red 
and blue galaxies on the red sequence is an incomplete 
model. In order to match the variations of the spectral 
features across the red sequence, a more detailed model 
is required. This model should include both age and 
metallicity variation and, more importantly, non-solar 
abundance ratios. Until spectral synthesis models with 
realistic a-enhancement prescriptions are available, more 
quantitative analysis of our observed difference spectrum 
is difficult as the interplay between age, metallicity, and 
a-to-iron ratio variations must be included for a full 
understanding of the physical differences between these 
galaxies. 

Another clear difference between the models and the 
data is the strong Na D feature near 5892 A. Here, the 
red composite spectrum has considerably more Na D ab- 
sorption than predicted by the stellar synthesis models. 
As the model spectra do not contain absorption due to 
the interstellar medium, this discrepancy may indicate 
that the red galaxies contain more neutral gas than their 
blue counterparts. 

7. CONCLUSIONS 

We present new multi-color observations of a large 
sample of galaxies in the range 0.1 < z < 0.4 gathered 
from the Sloan Digital Sky Survey. These galaxies rep- 
resent the most massive galaxies in the Universe and are 
not selected to be in clustered environments. Past studies 
of the color-magnitude relationship have focused on L* 
galaxies in clusters, thus our work focuses on an area of 
parameter space previously unexplored. Using the robust 
photometry of bright galaxies provided by the SDSS as 
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Fig. 8. — (a) Composite spectra of galaxies between 1-a and S-a of the mean '^'^^{g — r) color. The gross properties of the composite red and blue 
spectra are quite similar, (b) The difference spectrum between the blue and red galaxies (middle) compared to the difference spectrum predicted 
for a passively evolving population (top) and an old population with a recent injection of young stars (bottom). The difference spectrum shows 
clear signs of enhanced Balmer absorption lines in the blue galaxies as well as stronger molecular features in the red galaxies. The red galaxies 
have stronger sodium absorption, leading to the "emission" feature in the difference spectrum, (c) Residuals between the data and comparison 
synthesis spectra. The general properties of the observed difference spectrum are matched in both comparisons, though the model consisting of an 
old galaxy with a small amount (~ 3%) of star formation in the past 1 Gyr slightly over-predicts the continuum strength in between 3900 A and 
4100 A compared to the passively evolving model. 



well as multiple observations of moderate-redshift galax- 
ies, we are able to construct the color-magnitude rela- 
tions for massive field galaxies in several colors. While 
our data only span one magnitude in luminosity, and thus 
are not ideal for in-depth analysis of the slope and zero- 
point of the CMR, our large sample of galaxies allows us 
to measure the dispersion about this relationship quite 
well. The scatter around the red sequence is 20 mmag 
in the reddest bands while the bluest colors have nearly 
double this value. 

We construct several toy model star formation histo- 
ries in order to place constraints on the evolutionary pro- 
cesses at work in massive galaxies. We find that the 
majority of star formation in these early- type galaxies 
must have occurred before z ~ 1, in agreement with 
other authors. For secondary bursts of star formation 
with strengths less than 5% of the total stellar mass of 
the system, our data are consistent with bursts of star 
formation at any epoch. However, one should note that 



none of the toy models reproduce the observed distribu- 
tion of galaxy colors. 

Furthermore, we make comparisons between our new 
measurement of the scatter around the color-magnitude 
relation for massive field galaxies and early-type galaxies 
in clusters. An internal comparison shows that the scat- 
ter around the CMR in dense environments is 11 ± 2% 
smaller than that measured using more isolated galaxies. 

The composite spectra of LRG galaxies redward and 
blueward of the mean galaxy color are physically dis- 
tinct. The blue galaxies show stronger Balmer absorp- 
tion than the red galaxies, consistent with the presence 
of a younger stellar population. Simple models slightly 
favor passive evolution of stellar populations as the cause 
of the observed differences in spectra across the red se- 
quence compared to a more recent injection of young 
stars onto an old stellar population, though the differ- 
ence between the predicted spectra from these models is 
small. 
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